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Abstract 

A new class of cosmological models in f(R, T) modified theories of gravity proposed by Harko et al. 
(2011), where the gravitational Lagrangian is given by an arbitrary function of Ricci scalar R and the trace 
of the stress-energy tensor T, have been investigated for a specific choice of f(R,T) — fi(R) + f2{T) by 
considering time dependent deceleration parameter. The concept of time dependent deceleration parameter 
(DP) with some proper assumptions yield the average scale factor a(t) = sinh~ (at), where n and a are 
positive constants. For < n < 1, this generates a class of accelerating models while for n > 1, the models of 
universe exhibit phase transition from early decelerating phase to present accelerating phase which is in good 
agreement with the results from recent astrophysical observations. Our intention is to reconstruct f(R,T) 
models inspired by this special law for the deceleration parameter in connection with the theories of modified 
gravity. In the present study we consider the cosmological constant A as a function of the trace of the stress 
energy-momentum-tensor, and dub such a model 'A(T) gravity" where we have specified a certain form 
of A(T). Such models may display better uniformity with the cosmological observations. The statefinder 
diagnostic pair {r, s} parameter has been embraced to characterize different phases of the universe. We also 
discuss the physical consequences of the derived models. 

Keywords: Bianchi type-V universe, Modified gravity, Variable deceleration parameter, Time dependent A-term 
PACS number. 04.50.Kd, 98.80.Es 

1 Introduction 

Despite the fact that Einstein's general relativity (EGR) is a triumphal success theory and the basis for the 
description of most of gravitational phenomena known to date, it fails to explain the recent discovery of the 
accelerating expansion of the universe. The recent observational data [1—13] strongly indicate that our universe 
is dominated by a component with negative pressure, dubbed as dark energy (DE), which constitutes with ~ 3/4 
of the critical density. In order to explain the nature of dark energy and the accelerated expansion, a variety 
of the theoretical models have been proposed in the literature such as quintessence, phantom energy, k-essence, 
tachyon, f-essence, chaplygin gas, etc. Among these different models of DE, the modified gravity theories are 
f(R) gravity [14—18] and Gauss-Bonnet gravity or f(G) gravity [19—23]. Another approach to modified gravity 
theory is so called f(T) gravity [24—26], where T is scalar torsion. These modified gravities have recently been 
verified to explain the late-time accelerated expansion of the Universe. 



One of the interesting and prospective versions of modified gravity theories is the f(R,T) gravity proposed 
by Harko et al. [27] , where the gravitational Lagrangian is given by an arbitrary function of the Ricci scalar R 
and of the trace of the stress-energy tensor T. They have obtained the gravitational field equations in the metric 
formalism, as well as, the equation of motion for test particles, which follow from the covariant divergence of 
the stress-energy tensor. The f(R,T) gravity models can explain the late time cosmic accelerated expansion of 
the Universe. Recently, Chaubey and Sukla [28 , Adhav [29], Samanta [3Q ; , Reddy et al. [31—33] have studied 
cosmological models in /(i?,T) gravity in different Bianchi type space-times. 
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In recent years Bianchi universes are playing important role in observational cosmology, since the WMAP 
data [34—36] seem to require an addition to the standard cosmological model with positive cosmological con- 
stant that bears a likeness to the Bianchi morphology [37—42]. According to this, the universe should reach to 
a slightly anisotropic special geometry in spite of the inflation, contrary to generic inflationary models [43—49] 
suggesting a non-trivial isotropization history of universe due to the presence of an anisotropic energy source. In 
order to explain the homogeneity and flatness of the presently observed Universe, it is usually assumed that this 
has undergone a period of exponential expansion [43, 45—47]. Mostly the expansion of the universe is described 
within the framework of the homogeneous and isotropic Friedman-Robertson- Walker (FRW) cosmology. The 
reasons for this are purely technical. The simplicity of the field equations and the existence of analytical solutions 
in most of the cases have justified this over simplification for the geometry of space-time. However, there are 
no compelling physical reasons to assume the former before the inflationary period. To drop the assumption 
of homogeneity would make the problem interactable, while the isotropy of the space is something that can be 
relaxed and leads to anisotropy. Several authors [50—54] have studied particular cases of anisotropic models and 
found that the scenario predicted by the FRW model stand essentially unchanged even when large anisotropies 
were present before the inflationary period. Thus, the anisotropic Bianchi models become of academical interest. 
Out of this, in the present study Bianchi type-V space-time is taken into consideration. 

Motivated by the above discussions, in this paper, we propose to study cosmological model in F(R, T) gravity 
in Bianchi type-F space-time by considering f(R,T) = fi(R) + f2{T). The outline of this paper is as follows: 
In Sect. 2, gravitational field equation in the f(R,T) modified gravity theory is established. Section 3 deals 
with the basic equations and quadrature solutions. In Sect. 4, the physical and geometric aspects of the model 
have been discussed. Section 5 deals with the physical acceptability of the solutions. Finally, conclusions are 
summarized in the last Sect. 6. 

2 f(R,T) modified gravity theory 

The theory suggests a modified gravity action given by 

S= ~ikJ f( R > T )^ d * x + J LmV=9d 4 x, (1) 

where f(R, T) is an arbitrary function of the Ricci scalar, R, and the trace T of the stress-energy tensor of the 
matter, T M „. L m is the matter Lagrangian density. The stress-energy tensor of matter is defined as 

2 S^gL 



Tfj, v — — . , (2) 

and its trace by T = g^ v T^ v . Assuming that the Lagrangian density L m of matter depends only on the metric 
tensor components g^ v and not on its derivatives leads to 

Tfu, = gfivL m - 2 q ■ (3) 

Here and in this subsection, we use the unit of G = c = 1. As past related studies, a theory whose Lagrangian 
density is described by an arbitrary function of R and the Lagrangian density of matter as F(R, Lm) has been 
explored in Harko and Lobo |55) . Moreover, in Poplawski |56j a theory in which the cosmological constant is 
written by a function of the trace of the stress-energy tensor as A(T) has been investigated. 

Varying the action S with respect to the metric tensor components g^ v , the gravitational field equations of 
f(R, T) gravity is obtained as 

f R (R, T)R^ - T)g„ v + (g MV n - V^V^MR, T) 

= 8ttT m1/ - f T (R, T)T^ V - f T (R, T)Q^. (4) 

with 6^ = gV (fpt), which follows from the relation 6 (j^) = T„ v + 9^ and □ = V'V,, f R (R,T) = 
eof^ ' It{R,T) = a ^^ ,T - ) and Vj denotes the covariant derivative. The contraction of Eq. Q yields 
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f R (R,T)R + 3Df R (R,T) - 2f(R,T) = (8tt - f T (R,T))T - f T (R,T)Q with 9 = g^Q^. Combining Eq. 
(jj) and the contracted equation and eliminating the nf R (R,T) term from these equations, we obtain 



f R {R,T)[R llv --Rg i 



:f(R,T)g^ 



i7T-f T (R,T))[T tlv --Tg liV 



f T (R,T)[e iil/ --Qg flu 



■V M V v /fl(/2,T). 



(5) 



On the other hand, the covariant divergence of Eq. ([T]) as well as the energy-momentum conservation law 
V M [f R (R,T) - \ f(R,T)g^ u + (g^U - V M V A1 )/ fl (i?, T)] = 0, which corresponds to the divergence of the left- 
hand side of Eq. ([!]) , we acquire the divergence of T M „ as 



fr(R,T) 



8Tr-f T (R,T) 
In addition, from T^ v = g^L M - 2 (§f$r) we have 



(T^ + e„„)V In f T (R, T) + ve 



(6) 



5gf* v ySgt* 



( S 9ij 1 I r r ., 



dg^dg 1 - 



Using the relation g^ii. 

by 



-gijdjcrSlZ with 8^1 = l^ir, which follows from g i7 g 
Q fiu — -)- g^i/Ljn 2^ 



7J 



(7) 



6j , we obtain M1 , as given 



(8) 



Provided that matter is regarded as a perfect fluid, the stress-energy tensor of the matter Lagrangian is given 
by 



(p + p)u^u„ -pg^ v , 



(9) 



where it M = (0, 0, 0, 1) is the four velocity in the moving coordinates which satisfies the conditions u^u v = 1 and 
u^WvUfi — 0. p and p are the energy density and pressure of the fluid respectively. With the use of equation 
flSJ, we obtain for the variation of the stress-energy of a perfect fluid the expression 



B — — IT — 



(10) 



Since the field equations in f(R, T) gravity also depend on the physical nature of the matter field (through the 
tensor M „), for each choice of / we can obtain several theoretical models. Three explicit specification of the 
functional form of / has been considered in Harko et al. [27] 

r R+2/(T) 

f(R,T)={ h(R) + f 2 (T) 

{ MR) + f2(R)f*(T) 

The cosmological consequences for the class f(R,T) = R + 2/(T) have been recently discussed in details by 
many authors [28—33, 57]. Recently, Shamir et al. [S5J and Chaubey & Sukla [55] have discussed Bianchi type-I 
& V and a general class of Bianchi models respectively in F(R, T) gravity by considering f(R, T) = R + 2/(T). 
In this paper we are considering the cosmological consequences of the class for which f(R,T) = fi(R) + fiiT), 
Our derived cosmological model is totally different and new from that of other authors mentioned here. So far 
the physically important cosmological term A which is a candidate for dark energy remains less attended. So, 
our derived model may lead to better understanding of the characteristic of Bianchi type-V models. 

The gravitational field equation @ becomes 

f'i{R)R P u - \h{R)g^ + (g^a - v M v„)/i(ii) = 
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8vrT M „ + f 2 (T)T^ + (f' 2 (T)p + -/aCo) g^ (11) 

where the prime denotes differentiation with respect to the argument, the field equations of the standard f(R) 
gravity can be recovered for p = (the dust case) and f 2 (T) = 0. We consider a particular form of the functions 
fi(R) = XiR and f2(T) — X 2 T where Ai and A2 are arbitrary parameters. In this article we take Ai = A2 = A 
so that f(R,T) = X(R + T). 

Equation (jlll) can now be rewritten as 

XR^ - ~X(R + T)g^ + (g^n - V M V„)A 

= 8ttT^ - XT^ + A(2T M „ + pg^). (12) 

Setting {g^a - V A1 V iy )A = 0, we get 

XG„ V = 8nT^ + XT^ + (x P + ^XT^j 9tiv , (13) 
where G M „ = i? M „ — ^g^R is the Einstein tensor. This could be rearranged as 

Gfiv -{p + \t\ g^ = ^^T^. (14) 

Recalling Einstein equations with cosmological constant 

G^y - Kg^ = SttT^. (15) 

We choose a negative small value for the arbitrary A so that we have the same sign of the RHS of (fT3"|) and (fl4"|) , 
we keep this choice of A throughout. The term (p + |T) can now be regarded as a cosmological constant. So, 
in the framework of the f(R,T) gravity, we can get the cosmological constant as a function of the equation of 
state parameter lo , the energy density p and the trace of the stress-energy tensor. But since lo and p are already 
included in T so we could just write 

A = A(T) = p + tjT. (16) 

The dependence of the cosmological constant A on the trace of the energy momentum tensor T has been proposed 
before by Poplawski [5 6) where the cosmological constant in the gravitational Lagrangian is a function of the 
trace of the energy-momentum tensor, and consequently the model was denoted "A(T) gravity" . It was argued 
that recent cosmological data favour a variable cosmological constant which is consistent with A(T) gravity 
without the need to specify an exact form of the function A(T). A(T) gravity is more general than the Palatini 
f(R) and could be reduced to it if the pressure of matter is neglected (Magnano|59j. Poplawski [60 ). Considering 
the perfect fluid case, the trace T = —3p + p, for our model. 

K= l -{ P -p). (17) 



3 The basic equations and quadrature solutions 

Modern observations (WMAP data for example) indicate that the universe is not completely symmetric [34—36]. 
From that point of view Bianchi models (which represents spatially homogeneous and anisotropic spaces) are 
more appropriate in describing the universe as it has less symmetry than the standard FRW models. We use 
the following metric of general class of Bianchi cosmological model 

ds 2 = dt 2 - A 2 dx 2 - e- 2 ^ x [B 2 dy 2 - C 2 dz 2 }, (18) 

where /3 is a constant and the functions A(t), B(t) and C(t) are the three anisotropic directions of expansion in 
normal three dimensional space. Those three functions are equal in FRW models due to the radial symmetry 
and so we have only one function a(t) there. The average scale factor a, the spatial volume V and the average 
Hubble's parameter H are defined as 

a = (ABC)^ (19) 
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V = a 3 = ABC, (20) 

and 

H=^(H X + H 2 +H 3 ), (21) 

respectively with Hi — ^ , H% — and H 3 = ^ . Here and elsewhere the dot denotes differentiation with respect 
to cosmic time t. From Eqs. (fT]), ^ and ^ we get 

IV 1 ( A B d\ 

h = zv = z{a + b + c)- (22) 

Now the cosmological equations for the energy momentum tensor ([9]) and the metric (|18[) are 

>-A, (23) 

p - A, (24) 
p - A. (25) 





C 


P 2 


B 


+ c~ 


A 2 ' 


A 


c 


P 2 


A 


+ c~ 


A 2 


A 


B 


P 2 


A 


+ B- 


A 2 ~ 


1 


BC 


3/3 2 


? + 


BC ~ 


A 2 




4- 


B 






B ~ 



A 



= 0. (27) 
In the next proceeding of the paper we take constant (/3) as unity without any loss of generality. 

Integrating Eq. (|2"7| and absorbing the integration constant into B or C, we obtain 

A 2 = BC, (28) 



without any loss of generality. Subtracting (123)) from (f24j) . (f23j) from (|25|) . and (j24|) from (j25|) and taking second 
integral of each, we obtain the following three relations respectively: 

^=4 |), (29) 

^ = d 2 exp ^fc 2 y ^ , (30) 

and 

^ = d 3 exp ^/c 3 y ^ , (31) 

where d\, g? 2 , d 3 , ki, fc 2 and k 3 are constants of integration. Finally, using a — (ABC)i , we write the metric 
functions from (|29 |) — (|3"Tj) in explicit form as 



A(t) = Ziaexp ^mi J a 3 dt^j , 
B(t) = Z 2 acxp ^m 2 J a~ 3 dt 



(32) 
(33) 

C(t) = Z 3 aexp \ m 3 J a~ 3 d?j , (34) 



where 



h = Vd&, h = \/d[ 1 d 3 , l 3 = i/(d 2 d 3 )-\ (35) 
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and 



ki + k 2 k 3 -ki -(k 2 + k 3 ) 

m i = o ; m i = n ) TO 3 = 



3 ' "" 3 ' ""° 3 
The constants 7711,7712,7713 and ii, ia s '3 satisfy the fallowing two relations: 

mi + m 2 + m 3 — 0, hhh = 1- 
Substituting Eq. ([28]) in Eqs. (|32]>- (|34|) . we obtain 

il = 1, Z2 = = Cl, 777i = 0, 7772 = -7773 = C2, 



(36) 
(37) 
(38) 



where c\ and C2 are constants. Again, substituting Eq. (|38p in Eqs. (|32|) — ([34 | . the quadrature form of the 
metric functions in terms of average scale factor a can be written as 



A(t) 



B(t) = c\a exp I c 2 



C(t) = — exp 

Cl 



dt 
dt 

C 2 I T 



(39) 
(40) 

(41) 



Following Pradhan 61 , Pradhan et al. [62] and Chawla et al. [63], the modified gravity held equations can 
be solved by considering the time-dependent deceleration parameter which yield the average scale factor as 



a(t) = (shah (at) ) > 



(42) 



The derivation and the motivation to choose such scale factor has already been described in details by Pradhan 
[BT] . This relation (|42|) generalizes the value of scale factor obtained by Pradhan et al. [64] [65] in connection 
with the study of dark energy models in Bianchi tyj>e-VIo space-time and cosmological models with variable q 
and A-term respectively. In literature it is common to use a constant deceleration parameter [28—33, 66—68], 
as it duly gives a power law for metric function or corresponding quantity. The motivation to choose such time 
dependent DP is behind the fact that the universe is accelerated expansion at present as observed in recent 
observations of Type la supernova (Riess et al. (4][5]; Perlmutter et al. pQ; Tonry et al. [9]; Clocchiatti et al. 
[10] ) and CMB anisotropies (Bennett et al. [B]; de Bernardis et al. [6U] : Hanany et al. [7D]) and decelerated 
expansion in the past. Also, the transition redshift from deceleration expansion to accelerated expansion is 
about 0.5. Now for a Universe which was decelerating in past and accelerating at the present time, the DP must 
show signature flipping (see the Refs. [71—73]. So, there is no scope for a constant DP at present epoch. So, in 
general, the DP is not a constant but time variable. 



Eqs. (f3"9"| — ([4"T]) with the help of (14"2"1) . we obtain the metric functions as 

A(t) = sinh»(crf)), 
"ca(-l)^ 



B(t) — ci sinh" (at) exp 



C(t) — — sinh" (at) exp 
ci 



2a 



c 2 (-ir 



2a 



h(at)F(t) 



■ cosh(at)F(t) 



where 



Fit) =1 + 1(1 + -) cosh 2 (crf) + i- ( 1 + £ ] [ 1 + I ) cosh 4 (at) + o[cosh(crf)] 6 
6\ 77/ 40 V 77/V 



(43) 
(44) 

(45) 
(46) 
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4 Physical and geometric properties of model 



Using Eqs. (|4*3"1) — (l4"5j) in ([25)1 — ([25)1 and solving with (IT71) . we obtain obtain the expressions for pressure (p), 
energy density (p) and cosmological constant (A) for the model (fT5|) as 

p{t) = — \- fd cosh 2 (at) + C 2 cosh 4 (at) + C 3 cosh 6 (at) + C 4 cosh 8 (at) 

-Tl it) V 



+ C 5 sinh H ^ ii (at) + C 6 ). (47) 



p{t) = -—^—[Kx cosh 2 (at) + K 2 cosh 4 (at) + if 3 cosh 6 (at) 



+ Ki cosh 8 (at) + K 5 sinh 11 ^ (at) + K 6 J . (48) 
(#1 - ACi) cosh 2 (at) + (K 5 - AC 5 ) sinh^ 11 ( at ) + (Jf 6 _ C* 4 )j , (49) 



A = 



2F 1 (t) 
where 

Fi(t) = 16n 2 (87r + A) sinh 2 (at), (50) 

d = An 2 c 2 2 + 12nc 2 , (51) 

C 2 = 18ncl + 3n 2 cl + 9cf , (52) 

C 3 = 2n 2 cj + 18c 2 ,, (53) 

C* 4 = -n 2 c^ + 6ncl + 9c\, (54) 

C 5 = IQn 2 , (55) 

C 6 = 4n 2 c 2 - 16a 2 n (56) 

K x = An 2 c 2 2 \ + Ylnc\\ - 96Aa 2 - 768aV (57) 

K 2 = 18nc^A + 3n 2 c^A + 9c^A, (58) 

K 3 = 2c 2 2 \n 2 + l&cjX, (59) 

K 4 = -n 2 c 2 A + 6nc^A + 9c 2 A, (60) 

K s = 16n 2 (32vr + 5A), (61) 

K 6 = 4n 2 \4 + 16n\a 2 + 356n7ra 2 . (62) 

Figure 1(a) depicts the variation of pressure versus time for A = —0.1, c\ = 1.1, c-i = 1, a = 1 and n = 0.1 as 
a representative case. From the figure we observe that pressure is increasing function of time. It starts from a 
large negative value and approaches to a small negative value near zero. From the discovery of the accelerated 
expansion of the universe, it is generally assumed that this cosmic acceleration is due to some kind of energy- 
matter with negative pressure known as 'dark energy'. 

The energy density has been graphed versus time in Fig. 1(b). It is evident that the energy density remains 
always positive and decreasing function of time and it converges to zero as t — > oo as expected. 

The behaviour of the universe in this model will be determined by the cosmological term A, this term has the 
same effect as a uniform mass density p e // = — A which is constant in time. A positive value of A corresponds 
to a negative effective mass density (repulsion). Hence, we expect that in the universe with a positive value 
of A the expansion will tend to accelerate whereas in the universe with negative value of A the expansion will 
slow down, stop and reverse. Figure 1(c) is the plot of cosmological term A versus time. From this figure, we 
observe that A is decreasing function of time t and it approaches a small positive value at late time (i.e. at 
present epoch). Recent cosmological observations [1—5, 9, 10] suggest the existence of a positive cosmological 
constant A with the magnitude A(Gh/c 3 ) « 10~ 123 . These observations on magnitude and red-shift of type 
la supernova suggest that our universe may be an accelerating one with induced cosmological density through 
the cosmological A-term. Thus, the nature of A in our derived models are supported by recent observations. 
For detailed description via different theoretical models and cosmography tests, one can see the recent review 
by Bamba et al. [74 where it is shown that dark energy universes may mimic the ACDM model currently 
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(consistent with the recent observational data). 



The physical parameters such as directional Hubble parameters (Hi), Hubble parameter (H), expansion 
scalar (9), spatial volume (V), anisotropy parameter (A m ), shear scalar (a) and deceleration parameter (q) are, 
respectively, given by 

a 

H x = - coth(ai), (63) 



n 



a Co 

H 2 = - coth(at) + 4- j (64) 



77 



sinh" (at) 



H 3 = - coth(crf) °4 , (65) 



77 



sinh » (at) 



9 = 3H = — coth(crf), (66) 
n 

V = ABC = (sinh(af ) ) » , (67) 




nc 2 



a coth(at) sinh ™ (at) _ 



(68) 



° 2 = < W 

ysmh" (at) ) 

q = n(l- tanh 2 (at)) - 1. (70) 

From the Eq. J7DJ , we observe that q > for t < ^ tanh _1 (l - i)i and q < for t > ^ tanh _1 (l - i)i. It 
is also observed that for for < n < 1, our model is in accelerating phase but for n > 1, our model is evolving 
from decelerating phase to accelerating phase. It follows that in our derived model, one can choose the value of 
n which gives the physical behavior of DP consistent with the observations. Also, recent observations from SNe 
la, expose that the present universe is accelerating and the value of DP lies to some place in the range —1 < q < 0. 

For the present Universe (t = 13.8 GYr) with q = -0.73 (Cunha [75]). Eq. ([70]) yields the following 
relationship between the constants n and a: 



■ tanh 1 



13.. 



0.27 



(71) 



It is self explanatory from the above relation that for the present Universe, the model is valid only for n > 0.27. 
It is observed that for < n < 1, our model is in accelerating phase but for 77 > 1, our model is evolving from 
decelerating phase to accelerating phase. It follows that in our derived model, one can choose the value of n 
which gives the physical behavior of DP consistent with the observations. 



From Eqs. (fBT)) and (JHHJ> , we observe that the spatial volume is zero at t = and the expansion scalar is 
infinite, which show that the universe starts evolving with zero volume at t — which is big bang scenario. 
From Eqs. ([43]) — (|45[> . we observe that the spatial scale factors are zero at the initial epoch t = and hence 
the model has a point type singularity [76]. We observe that proper volume increases exponentially with time. 
The physical quantities isotopic pressure (p), proper energy density (p), Hubble factor (H) and shear scalar (a) 
diverge at t = 0. As t — > 00, volume becomes infinite where as p, p, H, 9 approach to zero. It is important to 
note here that lim t ^o (35) spread out to be constant. Therefore, the model of the universe goes up homogeneity 
and matter is dynamically negligible near the origin. This is in good agreement with the result already given by 
Collins [77]. Finally, we can say that the model represents a shearing, non-rotating, expanding and accelerating 
universe, which starts with a big bang singularity and approaches to isotropy at the present epoch. 

The dynamics of the mean anisotropic parameter depends on two constant n and c 2 ■ From Eq. (|68l) . we 
observe that at late time when t — > 00, A m — > 0. Thus, our model has transition from initial anisotropy to 
isotropy at present epoch which is in good harmony with current observations. Figure 1(e) depicts the variation 
of anisotropic parameter (A m ) versus cosmic time t. From the figure, we observe that A m decreases with time 
and tends to zero as t — > 00. Thus, the observed isotropy of the universe can be achieved in our model at 
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present epoch. In the frame of pure f(R) gravity, the pioneer works (Nojiri and Odintsov [75J EHI) present very 
natural unification of the early-time inflation and late-time acceleration. According to this observation which a 
is sub-case of current model there is no place for large anisotropy after inflation. But the case in the present 
model is different. In our derived model q < for t > t c , where the critical time t c — itanh~ 4 (l — — )s. Also 
we observe that the average scale factor a(t) — > as t —> and also a(t) —> oo as t —> oo. This indicates that 
there exists inflation. 



The Ricci scalar for the solution is given by 

1 



R 



4 [cj(3n + n 2 ) + 24a 2 ] cosh 2 (at) 



8n 2 [cosh(at) - 1] 
+3cl(n 2 + 9) cosh 4 (at) + 2c 2 (n 2 + 9) cosh 6 (at) 
~ 2 (n 2 - 6n - 9) cosh 8 (at) + 48n 2 sinh^(at) - 48a 2 n + 4n 2 c 2 



(72) 



The evolution of Ricci scalar R with cosmic time is shown in Fig. 1(d). We observe that the curvature is positive 
through the whole evolution of the universe and R — ¥ as t — ¥ oo and R oo when t — > indicating the initial 
singularity. 



5 Physical acceptability of the solutions 



For the stability of corresponding solutions, we should check that our model is physically acceptable. 
Sound speed: 

It is required that the velocity of sound v s should be less than velocity of light (c). As we are working in the 
gravitational units with unit speed of light, i.e. the velocity of sound exists within the range < v s — (^^j ^ !• 



We obtain the sound speed as 



where 



and 



dp = £(*) 
dp h(t) ' 



g(t) = Asinh" (at) 



3c 2 2 (n + 3) 2 cosh 8 (at) + 24c^n + 3) cosh 6 (at) 
+9C2 (n 2 + In + 7) cosh 4 (at) + 6c 2 (n 2 + 6n + 3) cosh 2 (at) 

16n sinh 2 (at). 



12nc2 — 16na 2 



h(t) = sinh" (at) Qc\\(n 2 + 3) cosh 2 (at) + 18c 2 A cosh 2 (at) 
+9c?,A(n 2 +9)cosh 4 (at) + (n + 3) cosh 6 (at) + 3c?,(n + 3) 2 A cosh 8 (at) 
+ 256a 2 7r(n - 3) + 16Aa 2 (n - 6) + 512mr + 80nA 



(73) 



(74) 



(75) 



Here we observe that v s < 1. Fig. 2(d) depicts the plot of sound speed with time. We observe that v s < 1 
throughout the evolution of the universe. 
Energy conditions: 

The weak energy conditions (WEC) and dominant energy conditions (DEC) are given by 
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Figure 2: The plots of energy conditions, sound speed v s verses time and also s versus r. Here A 
Ci = 1.1, c% = 1, a = 1 and n = 0.1. 
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(i) p>0, (ii) p-p>0 and (iii) p + p > 0. 

The strong energy conditions (SEC) are given by p + 3p > 0. 

The left hand side of energy conditions have been graphed in Figs. 1(6) & 2 (a, b, c). From these figures, we 
observe that 

• The WEC and DEC for the derived model are satisfied. 

• The SEC is also satisfied. 

It has been shown by Wald [80] that under very general conditions all Bianchi cosmologies (except Bianchi 
IX) with a cosmological constant and an energy momentum tensor satisfying the strong and dominant energy 
conditions, will unavoidably enter a phase of exponential expansion. With the help of this result Jensen and 
Stein-Schabes [81] showed that if the number e-folds the Universe expands during its exponential phase is given 
by N then it will take a time of the order t ~ e 2Ar v / A for anisotropy to have any effect on the observable universe. 
One remarkable result is the independence of this result from the type or magnitude of the initial anisotropy. 
Statefinder diagnostic: 

The various cosmological observational data support the A cold dark matter (ACDM) model, in which the 
cosmological constant A plays a role of dark energy in general relativity. At the current stage, the ACDM model 
is considered to be a standard cosmological model. Sahni et al. |82j have introduced a pair of parameters {r, s}, 
called Statefinder parameters. In fact, trajectories in the {r, s} plane corresponding to different cosmological 
models demonstrate qualitatively different behaviour. The statefinder parameters can effectively differentiate 
between different form of dark energy and provide simple diagnosis regarding whether a particular model fits 
into the basic observational data. The above statefinder diagnostic pair has the following form: 

to differentiate among different form of dark energy. Here H is the Hubble parameter and q is the deceleration 
parameter. The two parameters {r, s} are dimensionless and are geometrical since they are derived from the 
cosmic scale factor a(t) alone, though one can reproduce them in terms of the parameters of dark energy and dark 
matter. This pair provides information about dark energy in a model— independent way, that is, it categorizes 
dark energy in the context of back-ground geometry only which is not dependent on theory of gravity. Hence, 
geometrical variables are universal. 

For our model, the parameters {r, s} can be explicitly written in terms of t as 

1 



[nH 2 ~ 3a 2 csch 2 (at) + a 3 coth 2 {l - 2csch 2 (at)}] , (77) 
_2 

s = — TTTi 7 \3a 2 csch 2 (at) - a 3 coth 2 {l + 2csch 2 (at)}] . (78) 

6nH z (2q — 1) 

So we can get the relation between r and s as 

s = —^-(r-l). (79) 
3(3-5) 

Figure 2(e) depicts the variation of s against r. From Fig. 2(e), we observe that s is negative when r > 1. The 
figure shows that the universe starts from an Einstein static era (r — > 00, s — > —00) and goes to the ACDM 
model (r = 1, s = 0). 

Therefore, on the basis of above discussions and analysis, our corresponding solutions are physically accept- 
able. 
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6 Discussions 



Evolution of Bianchi type-V cosmological model is studied in presence of perfect fluid and variable cosmological 
constant in f(R,T) theory of gravity (Harko et al. 2011), where the gravitational Lagrangian is given by an 
arbitrary function of the Ricci scalar R and of the trace of the stress-energy tensor T. For each choice of the 
function f(R, T), we get different theoretical models. Three examples have been given in Herko et al. (2011) (i) 
f(R,T) =R + 2/(T), (ii) f(R,T) = MR) + f 2 (T) and (hi) f(R,T) = h(R) + / 2 (i?)/ 3 (T). The cosmological 
consequences of the case (i) has been discussed before by several authors. In this paper, the gravitational field 
equation has been established by taking case (ii) into consideration. To find the deterministic solution, we have 
considered a time dependent deceleration parameter which yields a scale factor as a(t) = sirih** (at), where n 
and a are positive constants. For n > 1, this generates a transition of the universe from the early decelerating 
phase to the recent accelerating phase. 

The current observations of the large-scale cosmic microwave background suggest that our physical universe 
is expanding isotropic and homogeneous models with a positive cosmological constant. The analysis of CMB 
fluctuations may confirm this picture. But other analysis reveal some inconsistency. Analysis of WMAP data 
sets shows that the universe could have a preferred directions (Pullen et al. [83]; Samal et al. [84]; Groene- 
boom and Eriksen 85!; and Armendariz-Picdn et al. [5(5]). The ILC— WMAP data maps show seven axes well 
aligned with one another and the direction Virgo. For this reason Bianchi models are important in the study of 
anisotropies. 

The main features of the models are as follows: 

• In summary, we considered the modified gravity which naturally unifies two expansion phases of the 
universe: inflation at early times and cosmic acceleration at current epoch. 

• The models are based on exact solutions of the f(R, T) gravity field equations for the anisotropic Bianchi- V 
space-time filled with perfect fluid with time dependent A-term. 

• The model represents an expanding, shearing, non-rotating and accelerating universe. 

• For suitable choice of constants the anisotropic parameter A m tends to zero for sufficiently large time (Fig. 
1(e)). Hence, the present model is isotropic at late time which is consistent to the current observations. 

• It has been found that A is a decreasing function of time and it converges to a small positive value at late 
time. The nature of decaying vacuum energy density A(t) in our derived model is supported by recent 
cosmological observations. These observations on magnitude and red-shift of type la supernova suggest 
that our universe may be an accelerating one with induced cosmological density through the cosmological 
A-term. 

• In literature it is a plebeian practice to consider constant deceleration parameter. Now for a Universe 
which was decelerated in past and accelerating at present epoch, the DP must show signature flipping as 
already discussed. Therefore, our consideration of DP to be variable is physically justified. Our derived 
model is accelerating at present epoch. 

• For different choice of n, we can generate a class of viable cosmological models of the universe in Bianchi 
type-V space-time. For example, if we set n = 1 in Eq. ([39]) . we find a = smh(at) which is used by 
Pradhan et al. |64) in studying the accelerating dark energy models in Bianchi type- Wo space-time and 
Pradhan et al. 65J in studying Bianchi type-I cosmological models with time dependent A-term. It is 
observed that such models are also in good harmony with current observations. 

• We observe that our derived solutions are physically acceptable in concordance with the fulfillment of 
WEC, DEC and SEC. The sound speed v s < 1 (see in Fig 2). The models in f(R, T) theory of gravity has 
stability and has initial singularity. Wald [80 pointed out that Bianchi cosmologies (except Bianchi IX) 
with A and satisfying DEC and SEC will enter in a phase of exponential expansion. 

• {r, s} diagram (Fig. 2(e)) shows that the evolution of the universe starts from Einstein static era (r — > 
oo,s — > —oo) and approaches to ACDM model (r = l,s = 0). So, from the Statefinder parameter {r, s}, 
the behaviour of different stages of the evolution of the universe have been generated. 
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Thus, the solutions demonstrated in this paper may be useful for better understanding of the characteristic 
of Bianchi type-V cosmological models in the evolution of the universe within the framework of f(R, T) gravity 
theory. 
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